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Abstract 

A new class of explicit finite- difference schemes for the computa- 
tion of weak solutions of nonlinear scalar conservation laws is presented 
and analyzed. These schemes are uniformly second-order accurate and 
nonoscillatory in the sense that the number of extrema of the discrete 
solution is not increasing in time. 

1. Introduction 

It is well-known that TVD (total variation diminishing) schemes, 
see e.g., [4] and the references given there, can give high resolution to 
discontinuities of the solution. TVD schemes, however, have only first- 
order accuracy at extrema. Of all second-order accurate TVD schemes, 
the Minmod scheme is the most diffusive. In [2], Harten and Osher devise 
the concept of nonoscillatory piecewise-linear reconstruction and modify 
the Minmod scheme to the UN02 scheme. Although this latter scheme is 
uniformly second-order accurate, it retains some of the diffusive character 
of the former. In this paper, we introduce a new class of schemes, named 
SONIC, that are uniformly second-order accurate and nonoscillatory. 

In section 2, we present TVD and UN02 schemes from a geometric 
framework which is rather different from the original presentations. Our 
description and analysis of these schemes naturally lead to the SONIC 
schemes later in the section. In section 3, the SONIC schemes are applied 
to nonlinear scalar conservation laws. It is then shown that the conjecture 
of Harten and Osher that the midpoint rule is nonoscillatory is indeed 
correct not only for the UN02 scheme but also for the SONIC schemes. 
In section 4, we present some computational results. We note in passing 
that R. Sanders [3] has recently developed a third-order accurate TVD 
scheme which degenerates to second-order accuracy near extrema. Our 
approach, however, is different and appears to be simpler. 
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2. Advection equation 

Consider the initial value problem for the advection equation with 
constant speed a > 0: 

du du 

m +a si = 0 ’ (21) 

u(x,0) = uq(x). 


The initial data uq(x) are assumed to be piecewise smooth functions that 
are either periodic or of compact support. Let u” be an approximation 
to the solution u at xj = jAx and the n-th time level. Equation (2.1) 
can be discretized in conservation form using the midpoint rule 



Uj — a\(u 


n+l/2 

>+ 1/2 


— U 


n+ l/2x 
>- 1/2 ) 


( 2 . 2 ) 


where A = At /Ax and an approximation to u at xj + 1 / 2 = 

(j + l/2)Ax and t n+1 / 2 = t n + 1/2A<. 

Let u ’} +] /2 be an approximation of u(xj+\/ 2 ,t n ). The value Uj+lj.J 
can be approximated by assuming that u n is linear in \xj,Xj+\/ 2 \- (See, 
e.g., [5].) At time t n+1 / 2 , the profile propagates to the right a distance 
aAt/2. The corresponding face value is 


U j+l/2 — U ]+ 1/2 — r ( u >+l/2 — u j )• (2-3) 

where r = a(At/Ax). Using (2.3), equation (2.2) can be written as 
Uj +1 = Uj - T(uj+i/2 - Uj-i/l) + T 2 (Uj+i/2 - Uj - Uj-i /2 +Uj-i) (2.4) 


(The superscript n has been omitted here and will be omitted in the 
remainder of the paper.) Throughout this section, the CFL niunbcr is 
assumed to satisfy r < 1. Our main task is to define the face values 
Uj+ 1/2 from the data Uj. The solution can then be obtained from (2.4). 

The following observation is useful for later proofs: let x(r) = Xj — 
t A x and y(r) =RHS of (2.4). The points (x = .r(r),?/ = y(r)) lie on 
a parabola Pj as r varies (At varies, a and A.t are fixed). Clearly, 
the points Uj = ( Xj,iij ) = (r(0), y(0)) and Uj - 1 = (xj-i,u.j-\) = 
(,r(l), y(l)) are on this parabola. The tangent T } to Pj at Uj has slope 
(uj + 1/2 — u j—i / 2 )/ Ax, i.e., Tj is parallel to the line — 1 / 2 ^>- 1- 1 /2 where 

Uj- 1/2 = (xj-i/ 2 , Thus Pj is determined by Uj,Uj-i, and 

Tj. The solution u” +1 in (2.2) is simply the y — component of the point 
(x(r),?/(r)) on the parabola Pj where 0 < r < 1 as shown in Fig. 1. 


2.1 TVD schemes 

From a geometric point of view, we describe and prove some im- 
portant properties of TVD schemes. Let I(z \ , . . . , z*) be the smallest 
closed interval containing Observe that median(x, y, c) lies 


2 


in the interval defined by any two of the three arguments. Similarly, 
median(zi,Z2,Z3,Z4,Z5) lies in the interval defined by any three of the 
five arguments, e.g., 

median (x,y, z) € I(y, z)] median(zi, . . . , zs) € I(zi, 24,23) (2.5a, b) 

To carry out the analysis, it is convenient to define the following central 
and upwind face values 

“i+1/2 = \( u j + u i+i)> ( 2 - 6 ) 

u j+l/2 ~ u j "b ( u i ~ u j— 1/2)1 ( 2 - 7 ) 

«“+ 1/2 = u i + (Uj - Uj- 1 ). ( 2 . 8 ) 

Notice that the Lax-Wendroff scheme corresponds to Uj+ 1/ 2 = Uy+1/2 
and the Warming- Beam scheme, to Uj+1/2 = u ]+i/2- 




The Minmod scheme, which is the simplest and the most diffusive of 
all second-order TVD schemes, can now be described in our framework. 
The face value for this scheme is 


Uj+ 1/2 = median (uj, ^ +1/2 , u“+ 1/2 )- ( 2 -9) 

For the sake of explaining the terminology, one can observe that 
minmod(x, y ) = median(a:, y, 0), 

median(x, y, z) — x + minmod(y — x, z — x) = y + minmod(x — y,z — y). 
The Minmod scheme posesses the following upwind monotonic property 

u] +1 € J(tti,«*i-i). (2.10) 

To show (2.10), observe that from definition (2.9) and property (2.5a), 
Uj- 1/2 G I(uj~ i,Uj_i/ 2 ); Uj + 1/2 € 1/2)- 
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Since u C j_ t j 2 lies between uj- 1 and Uj, and j^ 2 lies between Uj and 

up 2 

U j+ 1/2’ 

Uj— 1/2 G I(uj—\,Uj)\ Ujj r \j 2 G I{uj,Uj_ £ 1 / 2 )- (2.11) 

Properties ( 2 . 11 ) lead to the upwind monotonic property ( 2 . 10 ) by the fol- 
lowing reasoning. Another way to state ( 2 . 11 ) is that the point y 2 = 
(xj-x/2, Uj-1/2) lies on the line segment (xj-i/2,Uj-i)(xj-i/2,Uj) and 
U j+ 1/2 lies on (xy+1/2, Uj)Uj£ J /2 where ^“+1/ 2 = (*7+1/2, «“+ 2 i/ 2 )- Since 
the points (xy_i / 2 , Uj-i ), Uj and ^ u +j/ 2 are on same line, the slope of 
Uj~i/2Uj + i / 2 lies between 0 and s 2 where s 2 = slope([ 7 jd 7 “^j/ 2 ). Because 
the tangent T ; of the parabola Pj at Uj is parallel to t/y —1/2 ^>-1- 1 /2 by the 
observation following ( 2 . 4 ), slope(Ty) also lies between 0 and s 2 (see Fig. 
2(A)). Let Q 1, Q 2 be the parabolas though Uj, Uj-\ such that at Uj, Q 1 
has horizontal tangent and Q 2 has tangent of slope s 2 . It can be shown 
that Q2 has horizontal tangent at Uj-\ (one can show this by using a 
linear change of coordinates such that Uj-i — (0,0) and Uj = (1,1)). 
Q 1, Q2 form the boundary of the shaded region in Fig. 2(B). Since Pj 
also goes through Uj- 1, Uj and at Uj, its tangent lies between that of Q j , 
Q 2 , one can derive that the parabola Pj lies in the above shaded region. 
This completes the proof. 

Roe’s Superbee scheme can be interpreted in our framework as fol- 
lows: The face value for this scheme is 

Wj+1/2 = median(uy,t^ +1/2 ,wy + i,u“^ 1/2 ,u“^i /2 ). (2.12) 

The Superbee scheme is also upwind monotonic, that is, it satisfies prop- 
erty ( 2 . 10 ). To show this, observe that from ( 2 . 5 b), 

u j- 1/2 € /(uy-i, Wy_ 1/2 ,«>); uy+ 1/2 € I(uj,u U jl l/ 2 ,u U jf l/2 ). (2.13) 

As in the Minmod case, «y_i/ 2 lies between Uj-\ and iij, u u -^_ 1 j 2 lies 
between Uj and u j j 2 , thus (2.10) follows. 

To describe TVD schemes in general, let 

r = (uy-tiy_ 1/2 )/(uy +1/2 -uy); R = (uj-Uj-i )/(uy+t -Uj), (2.14a, 6) 

(see [ 4 ]). From definition ( 2 . 6 ), r — R. (However, this will not be the 
case later when we redefine Uy +1 / 2 .) The face value uy+1/2 is obtained 
by using a limiter function (f> 

Uj+ 1/2 = Uj + 4 >(r)(u C j +1/2 - Uj). (2.15) 
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The first-order upwind scheme corresponds to 4> = 0, the Lax-Wendroff 
scheme to <f> = 1, and the Warming-Beam scheme to <f> = r. 





Next, it is shown that the fact that <f> satisfies the TVD condition is 
equivalent to the fact that the face value in eq. (2.15) satisfies properties 

(2.11) . Indeed, the first half of (2.11) with index j replaced by index j + 1 
is Uj + i /2 G I(u i ,u i+ i).This is equivalent to the condition 0 < (f>(r ) < 2 
for all r. The second half of (2.11) is equivalent to <f>(r) G 7(0, 2r). Thus 

(2.11) is equivalent to the condition that <f> lies in the intersection of 
the above two regions, which is the shaded region and the negative r- 
axis of Fig. 3(A). This is the TVD region as defined by Sweby [4]. Since 
the Lax-Wendroff and the Warming-Beam schemes are both second-order 
accurate, to assure second-order accuracy away from extrema, the face 
value must lie between those of Lax-Wendroff and Warming-Beam, i.e., 
<f>(r) G 7(1, r). The intersection of this region and the TVD region is 
named the TVD2 region and is shown in Fig. 3(B). Observe that the 
upper boundary of the TVD2 region corresponds to the Superbee scheme, 
and the lower, to the Minmod scheme. 


2.2 UN02 schemes 

To resolve the problem of first-order accuracy at extrema of TVD 
schemes, Harten and Osher developed the UN02 scheme in [2]. From 
our point of view, this scheme is the same as the Minmod scheme except 
that the central face value is now obtained by a nonoscillatory 

quadratic interpolation. 

Let Qj be the parabola through Uj-i, Uj, and Uj+\. Since Qj and 
Qj+i have two points in common, they are either identical or have no 
other point in common. In the interval [xj,Xj + i], consider the three 
curves Qj, Qj+i and the line segment UjUj+ 1. Let Q>+i/2 be the curve 
which lies between the other two (Fig. 4). Let Uj +1 / 2 = (x J+ 1 / 2 , uC j+ 1 / 2 ) 
be the intersection of the two tangents to Qj+ 1/2 at Uj and Uj+ \ (notice 
the x location). An easy calculation yields the central face value 

w >+i/2 = ~2^ u i u j+ 1 ) — ~^Dj+ 1/2- (2.16) 


where 


Dj = Uj-i + Uj + i — 2 Uj; Dj+i /2 = minmod(Dj, Dj+i), (2.17a, b) 
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Having defined u< j+ i/2’ the face value of the UN02 scheme can readily 
be obtained from eqs. (2.7) and (2.9). 



(A) (B) 

FIGURE 4. FIGURE 5. 


The first essential property of this scheme is that the interpolation 
(2.17), (2.16) is nonoscillatory in the following sense: If is a strict 

local maximum (similarly for minimum), i.e., Uj < Uj +1 / 2 a *id Uj + 1 < 
u °j+ 1 / 2 > then the data (uj.) has a strict local maximum at Xj or £j+i, 
i.e., Uj - 1 < uj and Uj + 2 < u j+i • To prove this, observe that since 
Uj+ 1/2 is a strict local maximum, the parabola Qj+ 1/2 is concave down 
with the maximum value lying between Xj and x ]+1 . At least one of the 
two points Uj— 1, Uj- \-2 lies on this parabola. If both of them belong to 
Qj+ 1 / 2 , the above claim follows. If only one of them lies on Qj+i/ 2 , then 
the other point belongs to the parabola Q* +1 / 2 which lies above Qj+ 1/2 
in the interval (xj,Xj+i) by definition. Since two different parabolas can 
have at most two points in common, Q*j +\/2 lies below Qj+ 1/2 outside of 
[xj,Xj+i] as shown in Fig. 5(A). This completes the proof. The above 
property implies that the interpolation procedure does not increase the 
number of strict local extrema. 

The next property of the UN02 scheme is that if Uj +1 / 2 is a strict 
local maximum (similarly for minimum), then 

u j+l /2 — u j+i/2> u j+ 3/2 = u j+3/2- (2.18a, b) 

Indeed, let Vj-\ = (xj-i,Vj-i), Vj+2 = (xj+2,Vj+ 2 ) be on the parabola 
Qj+1/2 (see Fig. 5(B)). From the above proof, Uj-i < Vj-i and Uj + 2 < 
Vj+2- Let Vj-i/2 — ( x j-i/2i v j-i/2) t> e on th* 3 H ne Uj Uj +1 / 2 . Then 
Uj_ j/2 < Vj_ x / 2 . This implies (2.18a). As for (2.18b), it follows from 

C ^ up 

u i+ 1/2 - u j+ 1/2- 

Let RQj-1/2 be the range of the parabola Qj-1/2 in [xj-i,Xj], that 
is, the interval formed by the values Qj_i/ 2 (x) where x varies in [xj-i , Xj\. 
We can now show that the UN02 scheme posesses the following quadratic 
upwind monotonic property 

u" +1 e I(uj,u c j_i/ 2 ,uj-i). (2.19) 

In fact, we will show u” +1 6 RQj-1/2 which implies (2.19). Observe 
first that (2.13) holds. If u C j _ 1 / 2 is not a strict local extremum then 
RQj-i/2 = I{uj-i , Uj) and (2.19) follows exactly as in the Minmod case. 
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If u j-i /2 i s a strict local extremum, say a maximum, then eqs. (2.18a,b) 
imply that the parabola Pj is precisely Qj- i/ 2 . This completes the proof. 

A quadratic upwind monotonic scheme is nonoscillatory in the sense 
that it can not create new strict local extrema. This follows from the 
nonoscillatory property of the interpolation (2.16), (2.17). 

We next show a simple proof of the second-order accuracy of the 
UN02 scheme. Denoting h = Ax , one obtains by using Taylor series 
expansions, 

(«,'+,/ 2 - «>)/(*>/ 2) = u'j + 0(h 2 ), (2.20a) 

(u, - u‘_ 1/2 )/(h/ 2) = + 0(h 2 ), (2.20 b) 

(*#,/» - Uj)/(h/2 ) = (uj - uj_ l/2 )/(h/2) = «' + 0(/> 2 ), (2.20c) 

(«> - u^ 1/2 )/(h/2) = + 0(h 2 ). (2.20 d) 

where it' denotes du/dx. The last expression follows by writing (2.20c) 
with index j replaced by j — 1 and then expanding u' _ j and itj_i around 
Xj. Since the face value it J+ i/ 2 must lie between u C j + i/ 2 a nd Uj+i / 2 » 

(w 2 +i/2 — Uj-\/2)/h = it' + 0(h 2 ). (2.21) 


Thus the scheme is second-order accurate. 

We have presented the UN02 scheme as a modification of the Min- 
mod scheme by redefining «^ +1 / 2 using eq. (2.16) above. Will the same 
procedure yield a modification of the Superbee scheme? The answer is 
yes as will be shown in the next section. 

2.3 SONIC (second-order nonoscillatory interpolating conser- 
vative) schemes 

We combine the names MUSCLE [6] with UN02 and come up with 
\ONIC. Of all SONIC schemes, the UN02 scheme above is the most diffu- 
ve and the SONIC-B or Sonicbee scheme below is the most compressive, 
tis scheme can be defined by eqs. (2.17a,b), (2.16), (2. 8), (2. 12), (2.3) 
d (2.2). It possesses all the properties of the UN02 scheme that are 
\itioned above. To prove the second-order accuracy of the Sonicbee 
We, consider the following cases at face j + 1/2: 
yCl: Uj +i /2 is a strict local extremum, then u>+i /2 = W |+l/2 ; 
v2: Uj_ i/ 2 is a strict local extremum, then Uj +i / 2 = x / 2 ', 

W if neither of the above cases hold and if uj is a local extremum, 
in other words, ( uj — — Uj +i ) > 0, then uj + j/ 2 = Uj G 

: if none of the above cases hold, then we are in the regular case 
of either ascending order Uj-i < ^j_i/ 2 < u j < u j+i /2 ^ u j+i 
or descending order. The case RGl corresponds to r = (uj — 

— 1 /2 )/ ( W >-(-l/2 — U j) — 
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RG2: this is the regular case where r > 1. 

In any of the above cases, Uj+ i/ 2 G I{u C j+ l / 2 , ^+ 1 / 2 ) an d the above 
claim follows. 

The Sonicbee scheme is also quadratic upwind monotonic and the 
proof is exactly the same as that of the UN02 scheme. To describe the 
class of SONIC schemes, we need some more definitions. 

At each face j + 1/2, let the linear upwind monotonic interval Ij+ 1/2 

be the intersection of I(uj, Uj+i) and I(uj, uJ+i/ 2 ). Then condition (2.11) 
is equivalent to uj. |_i/ 2 G lj+ 1 / 2 , and, as shown in section 2.1, these are 
equivalent to the TVD condition in [4]. Similarly, let the quadratic up- 
wind monotonic interval Ij+ 1/2 be the intersection of I(uj,u C j +l / 2 ,Uj+i) 
and /(+,,«“£ 1 / 2 ,u]+ 1/2 ). Then condition (2.13) is equivalent to 

Uj+1/2 € l] + i/2- ( 2 . 22 ) 

Schemes satisfying (2.22) posess the quadratic monotonic property (2.19). 
The proof is similar to that of the UN02 scheme and is omitted. Observe 
that Ij+1/2 is equal to Ij+ 1/2 except in the cases EX1 and EX2 above 
where lj+ x j 2 strictly contains Ij+i/ 2 - Enlarging Ij+i /2 to Ij+ 1/2 these 
cases permits second-order accuracy at extrema. 

Let Ij+ 1/2 — I( c ii c 2 ). Any high-order accurate scheme can be made 
nonoscillatory by defining 


«j+i /2 = median(u J+ i/ 2 ,ci,c 2 ) 


(2.23) 


where Uj +\/ 2 is the face value of the scheme. 

We can now describe the SONIC schemes. Consider first the EX 
cases. In each of these cases, both the most diffusive (UN02) and the 
most compressive (Sonicbee) members give the same face value. Thus 
all other SONIC schemes must also have the same face value. One can 
also come to the above conclusion by observing that the intersection of 
^( u j+i/2’ m J+i/ 2) ar *d Ij+ 1/2? he., the intersection of the second-order ac- 
curate interval and the quadratic upwind monotonic interval, contains ex- 
actly one point. This point must be the face value of all SONIC schemes. 

Next, consider the regular case. In the cases RGl, RG2 above, 

Ij+ 1/2 = I( u ji u j+ 1/2)’ Ij+ 1/2 = I( u ji u j+ 1) (2.24a, b) 

respectively. Let cf>(r) be a continuous function for r > 0. With uj +x / 2 
defined by eq. (2.15), using (2.23) and (2.24a,b), one obtains 

Uj+ 1/2 = media n(«j +1/2 ,+,, u“+ 2 1/2 ), (2.25a) 

+,+ 1/2 = median(u^ +1/2 ,Uj,u J+1 ) (2.25 b) 
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respectively. Schemes defined as above axe nonoscillatory. They are 
second-order accurate if their face values lie between u c j +l / 2 and u j+i/2i 
and this is the case if <f> is chosen so that <f>(r) E 7(1,7’). Observe that 
the UN02 scheme corresponds to <f>(r) = min(l,r) and the Sonicbee, to 
< i>(r ) = max(l,r). 

The SONIC scheme corresponding to vanLeer’s MUSCLE scheme 
can be obtained by simply defining <j>(r) = (1/2)(1 + r), or equivalently 

u< j+ 1/2 = 2^ u i+i/ 2 “7+1/2)- (2.26) 

We name it the SONIC- A scheme (“A” for “average”). 

We conclude this section with the following remark. In defining 
U /+l/2> the minmod function of eq.(2.17b) can be replaced by any TVD2 
limiter function. One then comes up with a class of schemes which 
posesses all the properties of SONIC schemes and admits SONIC schemes 
as a subclass. We will report on this enlarged class elsewhere. 


3. Conservation laws 

Consider the following scalar conservation law 


du df n 

1- — = 0 

dt dx 


(3.1) 


where / is a smooth convex function of u, that is, / lies above all its 
tangents. This implies that a(u ) = df/du is an increasing function of 
u. Let u be the unique sonic point of /, i.e., a(u ) = 0. Using the 
simple change of variables v = u — u and g{v) = f(u) — f(u), one may 
assume without loss of generality that u = 0 and /( 0) = 0. The simplest 
conservation law with all the above properties is the inviscid Burgers’ 
equation 


du 1 du 2 

dt + 2 dx 


(3.2) 


Equation (3.1) can be discretized in conservation form using the midpoint 
rule 


u”+> = u" - A (f H „ 2 - (3.3) 


where, again, A = A t/Ax and fj+ 1/2 = /(«>+ 1/2) and the tilde replaces 
the superscript n + 1/2, e.g., 1/2 = \'j\. For simplicity of the 

proofs, the CFL number \i is assxamed to satisfy 


sup(A|a(?i)|) = n < 1/2 


The main problem is to define the face value u J+1 / 2 . To avoid mix- 
ing up the numerics and the physics, we adopt the projection-evolution 
approach advocated by vanLeer [5]. 
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The projection stage is purely numerics. First, assuming that the 
flow comes from the left, using the SONIC schemes of the previous section, 
one can define face value uj +1 ^ 2 . Then, on reversing the flow direction, 
a is now negative, one obtains tzj* +1 / 2 . This is the end of the projection 
stage. 

In the evolution stage, physics enters. We have to decide the flow 
direction and define the face value Uj+1/2 properly. To do this, consider 
the following cases, (cf. section 6 of [2] and section 4 of [1]) 

1- af+1/2 > 0 where aj +1/ 2 = a(uf +1/2 ). This is equivalent to uf +1/ 2 > 

0. We have two subcases: 

a- a f +\/ 2 > 0, then upwind direction corresponds to j ; 

b. flj+1/2 < 0, we have a shock. 

2. < 0, equivalently uj ' +1 / 2 ^ 0- Again, we have two subcases: 

a. aj+1/2 < 0, then upwind direction corresponds to j -f 1; 

b. Qj+i/2 — 0, we have a sonic expansion fan. 

In order to prove the nonoscillatory property later, the definitions of 
Uj+i/2 below are somewhat different from those in the literature. 

In case la above, assume that the approximations for u and a are 
linear in [xj, xj+i/ 2 \- Then, using a characteristic argument, 

Q h 

U j+l /2 - U j+l /2 ~ 1 T A J +1/2 ~tH u j + i/2 - «>)• (3.4) 

1 -t- A{a,j +1 / 2 - a,j) 

The ratio in the above expression is a more accurate approximation to 
Uj+i/2- F° r Burgers’ equation, Uj + i / 2 is of the same sign as Uj + 1 / 2 . In 
general to assure that the sign has not been changed, with e a small 
positive number (say, 10“ 9 ), we define 

«j+ 1/2 = max(e,u J i +1/2 ) (3.5) 

u j+ 1/2 = Uj+1/2 (3-6) 

Case 2a is similar to the above case with proper sign changes, 

‘‘f-H/2 = “j+1/2 + 1 + A( _ a r) 2 2 + aj + i) A K R + ./2 - ■)• (3-7) 

“j+1/2 = min(-e,u J R +1/2 ) (3.8) 

Uj+l/2 — ^j+1/2 (3.9) 

In case lb, we define uf +1 / 2 and Uj t +1 / 2 as in eqs. (3.4), (3.5) and 
(3.7), (3.8) respectively. The upwind direction depends on the sign of 

A7+I/2 = (f(Uj + j/ 2 ) — f(Uj + \/2))/(Uj+l/2 ~ Uj+ 1/2) (3.10) 
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If -A; + 1/2 > 0, Uj +\/2 is obtained from eq. (3.G) and if Aj+ 1/2 < 0, from 
eq. (3.9). 

Case 2b for the sonic expansion fan is simple: 

«;+ 1/2 = 0. (3.11) 

The above definitions of face values are consistent in the sense that if 
Uj+ 1/2 > 0, then wj+j/2 > 0 and the flow comes from the left. Conversely, 
if Uj+i /2 < 0, then uf +l/2 < 0 and the flow comes from the right. Finally, 
Uj + 1/2 = 0 corresponds to a sonic expansion fan where Uj + X j 2 < 0 and 
u f+ 1/2 — 0. 

Next, we show that SONIC schemes defined as above are nonoscilla- 
tory. The case uj+ 1/2 — 1/2 is trivial. Assuming that Uj+ 1/2 7^ Uj- 1/2, 

let 

Aj = {f(uj+ 1/2) - f(uj-i/ 2 ))/(uj + i /2 - Uj- 1/2), (3.12) 

then eq. (3.3) can be written in a form similar to (2.2) 

«" +1 = «” - A>A(ti i+1/ 2 - Uj- 1 / 2 ). (3.13) 

If Aj = 0 there is nothing to prove. Consider the case Aj > 0. Similar 
arguments hold for Aj <0. To show that these schemes are nonoscilla- 
tory, it suffices to show that they have the quadratic upwind monotonic 
property (2.19) where upwind direction corresponds to index j — 1 since 
Aj > 0. Indeed, with the CFL condition < 1/2, we only have to show 

Uj- 1/2 € Wj+i/2 € I(uj,Ujl 1/2 ,uJl 2 1/2 ). (3.14a, b) 

Notice that (3.14) is the same as (2.13) except that Uj- 1/2 and Uj+ 1/2 
are replaced by uj- x / 2 and Uj + j/ 2 . We use the above restricted CFL 
condition here. In the nonlinear case, the solution is no longer given by 
the parabola Pj in the observation following (2.4) but by the straight line 
Uj- 1 / 2 U jj. \ / 2 • Condition (3.14a) follows from the definition of face values. 
As for condition (3.14b), if Uj+1/2 > 0, then the consistency arguments 
after eq. (3.11) assert that at face Xj+ 1/ 2 , the flow comes from the left, 
thus the above claim. Before considering the case Uj+i / 2 < 0, we make an 
essential observation. Condition (3.14b) can be replaced by the following 
condition and (2.19) still holds: 

Uj+1/2 e I(uj-i/2,Uj + i/ 2 )- (3.15) 

See Fig. 6(A). Next, if Uj+ 1/2 = 0, then from the consistency arguments, 
Uj + 1/ 2 < 0. In order that Aj > 0, Uj-1/2 > 0 (Fig. 6(B)). Since Uj+1/2 = 
0, it lies between u>-i/2 a nd uf +i / 2 which is (3.15). Finally, we consider 
the case Uj+1/2 < 0, i.e., in the control interval [xj_i/ 2 , .Tj + i/ 2 ] there is a 
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shock since iij-1/2 > 0 (Fig. 6(C)). The difficulty here is that Uj+1/2 = 
uf+i/ 2 and yet we must show that (3.15) holds true. We need another 
definition. 


UP2 



(A) (B) (C) 

FIGURE 6. 


For each u, let v be of opposite sign such that f(v) = f(u). In 
the case of Burgers’ equation, v = — u. Because of the assumption that 
Aj > 0, 

Vj-1/2 < fij+1/2 < 0 < Vj+1/2 < Uj_ 1/2. (3.16) 

Rewrite (3.12) and (3.13) with Uj+i/ 2 replaced by Vj+i/ 2 : 

Aj = (f (vj+1/2) ~ f(uj -1/2))/ (vj+1/2 — Uj-1/2), (3-17) 

u] +1 = «" - AjX(v j+ 1/2 - Uj-1/2 )• (3.18) 

Observe that Aj\ is within the above CFL bound since 0 < Aj < a.j-1/2. 
Thus instead of (3.15) it suffices to show that 

Vj+1/2 e I(uj-i/2,Uj +1 / 2 )- (3.19) 

If uj + 1/2 < 0, (3.19) follows from (3.16). If Uj' +1 / 2 > 0, we have a shock 
at face Xj+i/ 2 and since the flow comes from the right, 

~ ~ R ~ L ~ 

Uj+1/2 — Uj + 1/2 < 0 < Uj +1 /2 < Vj + i/2- 

The last inequality above and the last one in (3.16) implies (3.19). This 
completes the proof. 

4. Computational results 

In the first set of problems, we present computational results for ad- 
vection with constant speed. The initial profiles consist of a sin 2 wave, 
a square wave, a triangular wave and a semi-ellipse wave. Each wave 
contains 20 grid points on a uniform grid of 200 points. Using periodic 
boundary conditions, the profile is advected one period with CFL num- 
ber 0.5. The results after 400 time steps are shown in Figs. 7 (A) for 
Minmod, (B) for UN02, (C) for MUSCLE, (D) for SONIC-A, (E) for 
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Superbee, and (F) for Sonicbee schemes. The figures on the left corre- 
spond to TVD schemes, and those on the right, to SONIC schemes. One 
can observe that the SONIC schemes are more accurate than their TVD 
counterparts. The Minmod scheme, and to a lesser degree, the UN02 
scheme are somewhat diffusive. The Sonicbee scheme can resolve the 
square wave and the triangular wave very accurately, however, it shares 
the compressive character of the Superbee scheme in the semi-ellipse wave 
and it has an overshoot in the sin 2 wave. The SONIC-A scheme appears 
to be the best of the schemes tested. Although it does not resolve the 
square wave as sharply as the Superbee and the Sonicbee schemes, it 
gives accurate solutions to the other three waves. 





(E) SUPERBEE SCHEME. (F) S0N1C-B SCHEME. 

FIGURE 7. - COMPARISON OF TVD (LEFT) AND SONIC (RIGHT) SCHEMES. 

In the second set of problems, results for Burgers’ equation are pre- 
sented. The initial condition is a sine wave of 20 grid points in the interval 
[—1,1]. The dots represent the numerical solutions and the curves, the 
exact solutions. With CFL number 0.5 and periodic boundary conditions, 
highly accurate results at time t = 0.4 are obtained in Figs. 8 (A) for 
UN02, (B) for SONIC-A, and (C) for Sonicbee schemes. Observe that 
a shock has formed at a; = ±1 and the solutions have no oscillations. 
The SONIC-A scheme is slightly more accurate than the UN02 and the 
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Sonicbee schemes in this case. 



(A) UN02 SCHEME. (B) SONIC- A SCHEME. (C) SONIC-B SCHEME. 

FIGURE 8. - NUMERICAL SOLUTIONS FOR BURGERS' EQUATION, t = O.M. 


5. Conclusion 

We have introduced the concept of upwind monotonicity which was 
used to show that the new class of schemes, named SONIC, is nonoscil- 
latory and second-order accurate. 
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